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1. INTRODUCTION 

During the last decade, researchers associated with tha aarth rasourcas 
program of NASA hava baan working on tha p rob 1 m of crop acraaga ami pro- 
duction estimation using LAI® SAT data. LAI® SAT Is a naar-aarth orbiting 
satalllta equipped with a multi spectral scanner (IKS) which measures tha 
reflectance of a target scene In various wavelength bands. The measure- 
ment unit 1$ a 1.1 acre square plot of land called a pixel. In estimating 
acreages devoted to a specific crop of Interest In a scene, each pixel Is 
assigned either to the crop or to the class of other ground categories; the 
pixel classification Is based on Its spectral response, say a kxl vector 
measurement [Heydom. et al . . 1979]. 

A scene Image In the form of a false color composite picture 1$ con- 
structed using Its MSS data. Image analysis and pattern recognition techniques 
are used to correlate the spectral characteristics In the scene to the features 
on ground. An area segment of several square miles is generally required for 
an analyst to interpret Its scene Image and to delineate discernible patterns 
for identifying possible land-use and land-cover classes. The MSS data for 
pixels of a class are modeled by a multivariate distribution function. Dis- 
criminant analysis techniques are applied to classify the MSS data and to 
estimate, say, a crop acreage proportion in a segment (Odell, 1976). 

As usual, a number of pixels are sampled to estimate the distribution 
parameters for the distinct classes of pixels and to specify the classi- 
fication procedure. Sampled pixels are first required to be Identified 
and labeled by their classes on ground. Lack of adequate spectral dis- 
crimination between the classes, among others, may cause mislabeling of 
of some pixels, thus resulting In a biased estimate of the classification 
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parameters. 

Another sou ret of trror In tstl Mating • crop terttgt proportion Is 
tilt prtstnet of mixed plxtls In « scant. A pixel Is defined mixed If 
It Is a boundary unit consisting of artas fro* tort than oat category of 
land use. Otherwise, It Is to be called a pure pixel. Often no distinc- 
tion Is made In the handling of nlx^d and pure pixels In clustering and 
classification of MSS data for estimating a crop acreage proportion In 
a segment. Previous empirical studies conducted at the Johnson Space 
Center have shown that the treatment of nixed pixels as If they are pure 
causes an additional bias In crop acreage estimation (Carnes and Baird, 
1980). A large-scale application of LANDSAT data for wheat estimation In 
U.S. and U.S.S.R. Is described In the Proceedings of Technical Sessions, 

The LAC IE Symposium, NASA (1979). 

The problem of estimating the relative acreage of a specific crop In 
an area segment can be formulated as follows: suppose C} and Co denote the 
classes of pixels for the crop of Interest and the group of other ground 
categories, respectively, and Cm denotes the class of mixed pixels In the 
segment. Considering the segment size to be large, let t m be the proba- 
bility of a random pixel to be from C m and for a spectral measurement X, 
let pi be the conditional probability that It belongs to Ci , 1-0, 1, 
given that the pixel Is from either Ci or Cq (l.e., It Is a pure pixel). 
Suppose p M i 1$ the proportion of acreages in C m that art devoted to the 
crop of Interest. Thus If p Is the actual acreage proportion for the crop 
of Interest In the segment, then 


P ■ »mPml + Pi 


( 1 . 1 ) 
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Suppose p Is estimated by the relative frequency of the npw t pixels 
that art classified Into Ci using a sample-based classification procedure. 
It Is assumed that only pure pixels are sampled for training the classifier 
and these are subject to mislabeling as Ci or C<). Let Rj and Rq be the 
classification regions for C\ and Co, respectively. Define the random 
variable 


I(*) 


1 if XeRi 
0 If XcRq. 


Then the estimate of p Is given by 

where N Is the total number of pixels In the segment (N Is assumed to be 
large). If N m 1 s the number of mixed pixels, then (1.2) can be written 
as 


P • *mPml + U-»m)Pl 


(1.3) 


where 

*m " l’’ Pml * ^ $l<h) 

N-Nni 

h » U/(N-iW 3 E K?i ). (1.4) 

In this paper, we Investigate p for its bias and variance. In section 
2, we express these parameters in terms of the expected classification error 
rates, their variances and covariance, and the first two moments of p m i. 
Considering the linear discriminant function for the classification rule, 
the asymptotic first two moments of error rates and those of p m i are obtained 
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1n section 3. Certain numerical results are given In section 4. It Is 
shown that the disparity In mislabeling rates of the two classes C} 
and Cq has a significant effect on the error rates as well as on the bias 
and variance of p. On the other hand, mainly the bias, and nrt the variance 
of p Is affected significantly as the relative size of mixed pixel class, 
w m varies. 

2. BIAS AND VARIANCE OF p 


Suppose 8} Is the error of classifying a pixel from Ci Into Co and 
00 Is for a pixel from Co Into Ci. Of course, there Is no error com- 
mitted In the classification of a mixed pixel. Let 0 m i be the probability 
of classifying a mixed pixel into C*. Suppose the classification rule Is 
determined on the basis of sample means and covariance matrices obtained 
from a sample of n pure pixels of which nj labeled as C} and nj> ■ n-nj 
labeled as Co. Let 7} and 7 q be the sample means, and §i and Sq be the sample 
covariance matrices for the two groups of labeled samples. Suppose 0}, 

§0 and Omi are the estimates of 0i, 8q, and 0 m i» respectively, 
given the sample observations. Then these estimates can be written In terms 
of the conditional probabilities as follows: 

01 * P[XeRol? eC l . Ti, Si, 1-0, 1] 

00 “ P[XeRl|XcC 0 , Ti, Si, 1-0, 1] 

®ml * p CXeRi|XeCm . Ji , §i , 1-0, 1] 


( 2 . 1 ) 
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Bits 

For the expected value of pj, we have 
£(0l) - ECEtfllli, §1. 1-0. 1]] 

» ECPfkRilJl. h» 1-0. IH 
Due to (2.1). we can write 

£(Pi) ■ Pi(i-£(8i)) ♦ po£(®o)« ( 2 « 2 ) 

So the bias of given by E(f>i)-Pi* Is 

B(pi) - -PiE(8i) ♦ PoE(®o)« ( 2 - 3 > 

As pointed out In the appendix, the expected acreage of devoted to 
the crop of Interest Is half of Its total size so that p m i-.5 and the bias cf 

Pml* 

B(Pml) ■ E(«.l) -5. (2*«) 

Accordingly* It follows from (1.3) that the bias of p, 

B(P) - w mB(Pml) ♦ 0-»m) B(Pi) 
where B(pi) and B(p m i) are given by (2.3) and (2.4). 

Variance 

For the variance: of pj, we can write 

Var(pj) • E[Var[pi|Ti , §i, 1-0,1]] 

♦ VarfECPilT,. S,. 1*0,13) 
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Slnce the tntlrt segment data art classified to obtain $i» tht conditional 
variance of ft, given sample data. Is zero. Thus, the first term on 
the right side Is zero, and 

Var(pi) - Var(PCX«fti|Ti, §t, 1-0,1]). 

Again, It follows from (2.1) that 

Var(ft) - pf Var(ft) ♦ pj Var(e 0 ) - ^pq O>v(e lv e 0 ). (2.6) 

Because of (1.3), we have the variance of p given by 

V.r(p) . 4 VirlPm) *■ ViPlPj) ♦ CovlPj P*1>(2.7) 

where Var(pj) is as given in (2.6), Var(p M i) Is simply the variance of $ m \ % 
and Cov(ft, p m j) obtained using the conditional argument. Is given by 

Cov(pi, p m i) ■ -Pi Cov($i, ©nil) ♦ PQ Cov(6 0 , 8 m i) (2.8) 


3. LINEAR JiSCRIMINANT ANALYSIS 


As considered by Heydorn, et a1.(1979), we assume that Co and C\ 
have multivariate normal distributions: X ~ N^yi, I) If Xcft , 

1-0,1. Without (OSS of generality, let 



, E ■ I 


(3.1) 


a 2 " («i - eoi' 5" 1 («i - i*o) 


where 
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Suppose Is the probability of mis labeling of • pixel from Cf, 1*0,1, 
for the sample labeling procedure, which Is generally manual using visual 
Interpretation of a scene Image and knowledge of crop characteristics In 
the area. The Image analyst who makes the labeling decision for sampled 
pixels uses their spectral Information plus hi* a-prlorl crop knowledge 
which Is fairly reliable at a somewhat larger level (e.g., crop field) 
but not at the pixel level. Obviously, labeling of a pixel Is partly 
dependent upon Its spectral response and as such It should be taken Into 
account . 

In another paper related to this topic, Chhlkara and McKeon (1983) 
have proposed an approach to modeling ml sal location. In general, and have 
discussed analytically the linear discriminant analysis In the presence of 
mlsal location In training samples. In the present context, their model 
(b) can be considered suitable for the mislabeling of pixels. This is 
a trusted model for which thresholds are determined for the two classes 
of pure pixels to assign class labels Cq and C} to sampled pixels. Since 
the crop Information Is utilized, these thresholds should be class specific 
and be functions of mislabeling rates. Given ag and the labeling 
procedure can be considered as follows: 

Suppose Xj Is the component of measurement vector X In the first 
dimension along wi.ich the class means are aligned. For XcCg, label 
the pixel as Cq If X| < -a/2 + Zl-2aQ and as Cl with probability 
.5, otherwise; and for XeCi, label the pixel as Cj If X} > a/2 + 
and as Cq with probability .5, otherwise, where Zi_2<»q and *2°^ 
are the (l-2a^)- and 2u - percentage points of the standard normal 
distribution. Under this rule, a pixel has at most fifty percent chance of 
mlsal location. Of course, one can consider other than half for the 
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maximum probability of ml sal location , say u, where (Ku< 1 ; 

but this would raqulra tha usa of Zi-Oq/u and *a : ' u for th# percentage 

points so that tha mislabeling ratas rtmaln as spaclflad. 

ChMkara and Me Kaon (1983) glva tha mlxtura distributions of tha two 
classas raprasantad In tha labeled training samplas and obtain tha asymp- 
totic distribution of tha sampla-basad boundary. Thair approach Is similar 
to that of Efron (1975) and can ba axtandad to obtain asymptotic first two 
moments of all three estimators, ej, 8 q and e n i. As discussed by 
Efron, the optimum boundary (l.e., tha case of known parameters) for tha 
linear discriminant rule Is a plana perpendicular to xj-axls and Inter- 
secting It at point t. whereas whan the sample size n Is large, tha 
sample-based boundary Is a plana Intersecting xi-axls at point t-kIt, 
with normal vector at an angle da from the xj-axl", where dt and da 
represent small deviations. If Dg. Oj and 0 m are the respective distances 
of the first component means of Cq, C} and C m from the optimum boundary, 
then their corresponding ! stances from the sample-based boundary are 
do ■ (0 q + dt) cos da 
di ■ (Dj - dx) cos da 

dm ■ (Dm + dt) cos da (3.2) 

Then, 8 q * ♦(- do) «nd 8j • e(-di). In Appendix we show that when 
A < 3.5, the distribution for can be approximated by normal , with 
Its mean zero and variance In the first dimension, of * 2/3 + a^/12. 

Though It Is not discussed. Its variance In any other dimension Is 1. 

Thus, the variance along the d^- direction 1$ of cos 2 da+ s1n 2 da 
or of + (1-oj 2 ) sin 2 da » oj, say. Accordingly, have 

$ml - ♦(- < W 0 d)* 
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By Taylor series expansion, Ignoring nlfher than second order dif- 
ferential term. It can be shown that 

6q ■ a(-0 0 ) - ♦(-t' 0 ) dt ♦ (1/2) D 0 g(-0 0 )C(dr) 2 ♦ (d«) 2 ] 

Sj - •(.0 1 ) ♦ ♦(-0 0 ) dr ♦ (1/2) 0 X ♦(-0 1 )C(dt) 2 ♦ (d«) 2 ] 

0*1 * ( dt/ ‘l> * i(V°l) 

♦ (-iy°l) C(dt/oj) 2 ♦ (de/ai) 2 ] (3.3) 

In (3.3), we have Dq ■ t ♦ A/2, Oj ■ -t * A/2 and ^ ■ t (e.g., 
refer to the figure given In Efron, 1975). how by a straight forward 
extension of the discussion and results of Chhlkera and McKeon (1983) the 
asymptotic first moments of So* ^1* ®m can obtained as follows: 

E(s 0 ) " •(-°o) + ii OoS(-Oo)C°t 2 ♦ O 

E(0 X ) * #(-0!) ♦ ^ D l ^(-D 1 )Co T 2 ♦ (k-1) ojl 

E(0 ml ) • ♦(-0*) ♦ ^ 0*s(-0*)[o T * 2 ♦ (k-1) o tt * 2 ] 

var(e 0 ) • /(s 2 {-D 0 )[o T 2 + ^ofto/ ♦ (ic-l)o tt 4 ]] 

Var(Sj) - /(♦ 2 (-D 1 )[o T 2 ♦ Jn0f[o T 4 ♦ (k-l)o u 4 ]] 

V*r(0 ml ) " /(♦ 2 (- D *)f°T* 2 + ^r(°*) 2 ^ T * 4 + (k-l)o tt * 4 ]] 

Cov(S 0 , 0j) - j0(-O o )0(-Oi)C-o T 2 ♦ (DQDj/En) 

(o t 4 ♦ (k-l)o u 4 )] 


- 11 . 


Co*(« 0 . > ml ) • d»(-D*)«(-0 0 ) [(« t 2 /«i) ♦ (D*Oo/aio{) 

(a , 4 ♦ (k-l)«. 4 )] 

cov(«i. i ml ) • C-(«, 2 /«i) ♦ <o*o 1 /a>«f) 

(o t 4 ♦ (k-1).. 4 )] (3.4) 


where 

Oq ■ t + A/2 , Oj » -t ♦ a/2 i 0* »• t/ji 

« T * ■ o T /tfi , ««* m 4 , J c l ( 3 « 5 ) 

with 2/3 ♦ a 2 /12, and o T 2 and aj as tht asymptotic variances 
for the discriminant isundary, which Are given In equations (3.10) and 
(3.11) of Ohlkara and McKeon (1983). 

4. NUMERICAL RESULTS 

In this section we Illustrate the olas and variance of p numerically 
by considering k-2, a- 2, p-.5, and n-100. First, In Table 1, we give 
values of t. o t 2 and o w 2 associated with the sample-based discriminant 
boundary when a-2, pj-.5, .3 and the mislabeling rates, apO and <*o«0, 

.1, .2, .3, .4. These values are taken from Table 1 In Chhlkara and McKeon 
(1983) corresponding to their model (b), and are used here to compute asymp 
totlc first two momentf of 6 q, 8j and 6 n | as described In (3.4). 

It Is seen that the discriminant boundary point, t shifts to the left 
as oq Increases. This Is expected due to disparity In mislabeling rates 
for the training samples disfavoring Cq which Is centered to the left on 
Xj-axls. The variance o r 2 Increases and oJZ decreases as a 0 Increases. 
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Table 1: Values of t, o t 2 and aj for the Sample-based Boundary (A«2) 



( 0 . 0 ) 

0 

.42 

1.000 

1.360 

2.000 

2.190 

(. 1 . 0 ) 

-.19 

.09 

1.136 

1.308 

1.068 

.845 

(. 2 . 0 ) 

-.40 

-.17 

1.541 

1.717 

.747 

.488 

(. 3 , 0 ) 

-.65 

-.44 

2.473 

2.542 

.644 

.387 

(. 4 . 0 ) 

- 1.00 

-.82 

5.373 

5.178 

.773 

.515 
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In Table 2, we show the results on mean e .d variance for different esti- 
mators as well as give the bias of pj and that of p w i when Pi«.5 and n-100. 
Based on the values of o T 2 and o w 2 given In Table 1, similar results 
can be easily computed for the case of Pi«.3. It Is seen from these 
results that E(0o) and Var(©o) Increase, whereas E(©i) and Var(0i) 
decrease as ao increases, but £(0^1) decreases and Var(&i) Increases. 

Again, tils can be expected because of a shift to the left In the boundary. 

§0 and §1 are negatively correlated, but has positive correlation 
with each of 0o and 0}. Interestingly, the variance of Is 

A 

affected only slightly, though It Is considerably higher than those of 0o 
and ©i when the disparity between the two mislabeling rates Is small. The 
absolute bias Increases for each of pi and p m i, and so are their 
variances and covariances. 

Next, we combine the two estimators pj and p m i by considering the 
proportion of mixed pixels, * m » 0, .1, .3, .4, .5, and compute the 
bias and variance of p (Table 3). Both the bias and the variance Increase 
as ag Increases. But there Is an Interaction with respect to change 
In w m ; the absolute bias and variance first decrease then Increase 
as w ra varies from 0 to .5. However, there Is only a slight change 
In the variance due to change In 


-14- 


Table 2: Means* Variances and Covariances of 8(j. and and 
Biases, Variances and Covariances of pi» p n i. 


(A-2 

» Pl*«5, 

n-100) 




Parameter 


Mislabeling Rates (ao*<»i) 




(0.0) 



(.3.0) 

(.4,0) 

E(«0) 

.162 

.212 

.276 

.365 

.500 

E(ei) 

.162 

.119 

.084 

.052 

.026 

E(Onil) 

.500 

.422 

.346 

.256 

.151 

Var(§o) 

.0006 

.0009 

.0017 

.0035 

.0085 

Var(8i) 

.0006 

.0004 

.0004 

.0003 

.0002 

*Var(8 m i) 

.0025 

.0017 

.0021 

.0026 

.0032 

Cov(8q,8x) 

.0006 

-.0006 

-.0008 

-.0009 

-.0012 

Cov (®0*%l) 

.0012 

.0013 

.0019 

.0030 

.0052 

Covteo.emi) 

.0012 

.0009 

.0009 

.0008 

.0007 

B (p 1 ) 

0 

.047 

.096 

.156 

.237 

B(Pml) 

0 

00 

r^. 

o 

• 

-.154 

-.244 

-.349 

Var(pi) 

.0006 

.0006 

.0009 

.0014 

.0024 

Cov(pi, p m i) 

0 

.0002 

.0005 

.0011 

.0023 


*Var(p m i) - Var(3 ml ) 
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Table 3: Bias and Variance of p (p».5, a«2» k»2, n»100) 


¥ ra 


Mislabeling Rates ( ag , «i) 




nrur 




(.4,0) 



(0 Bias 




0 

0 

.047 

.096 

.156 

.237 

.1 

0 

.035 

.071 

.116 

.178 

.2 

0 

.022 

.046 

.076 

.120 

.3 

0 

.010 

.021 

.036 

.061 

.4 

0 

-.003 

-.004 

-,004 

.003 

.5 

0 

-.016 

-.029 

-.044 

-.056 



(11) Variance 

xlO- 4 



0 

6 

6 

9 

14 

28 

.1 

5 

5 

8 

14 

27 

.2 

5 

5 

8 

14 

27 

.3 

5 

5 

8 

14 

26 

.4 

6 

5 

9 

14 

26 

.5 

8 

7 

10 

16 

27 



.16. 


5. CONCLUSION 

We have Investigated theoretically the error In crop acreage estimation 
using Landsat data and the current methodology of MSS data processing and 
linear discriminant analysis. Labeling of pixels by an Image analyst Is 
modeled and the effect of mislabeling rates on the bias and variance of the 
crop proportion estimate discussed. In past. Investigators have assumed a 
random model for ml sal location In training samples (Lachenbruch, 1966 and 
McLachlan, 1972), which Is not applicable here. Lachenbruch (1974) has dis- 
cussed two non-random models which are similar to our proposed model. He, 
however, studied these models only In the context of Fisher linear discri- 
minant function, assuming equal a -priori probabilities and evaluated Its 
performance using a simulation study. Presently no assumption of equal 
a -priori probabilities Is made and the numbers of pixels labeled as Co and 
C} are. In fact, treated as random, as one would expect. Only the total 
sample size n Is assumed fixed. 

This study extends the usual two-class classification methodology to 
a third class which presently arises due to mixed pixels In an area segment. 
Similar situation may also arise In inventorying forest, range, and other 
land-use and land-cover categories using a fallible measuring device. 

Presently we have assumed that the class of mixed pixels 1$ separable 
from the other two classes, and hence, N m Is known. If N m Is unknown and 
mixed pixels are delineated using an Imperfect boundary detection method, 
then an estimate of p Is obtained by 


where «m Is an estimate of tm. Because the number of mixed pixels In a 
segment Is large and because all such pixels are delineated (l.e., no sampl- 
ing Is envoi ved In the estimation of w m ), the variance of will be negllglbl 
Thus, we may consider * m known, assuming that the procedurt of delineating 
mixed pixels Is unbiased. 
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APPENDIX 

Distribution of Cm 


Let U be the proportion of acreage devoted to the crop of Interest In 
a randomly selected mixed pixel. Then U has the uniform distribution over 
Interval (0, 1). Thus, p ra j»E(U)».5. 

Dana (1982), and Lambeck and Potter (1979) have shown that the radiance 
received by Landsat sensor over a target area (resolution element) Is almost 
a linear function of the reflectivity directly transmitted from the target 
to the sensor. The aerosol optical thickness In atmosphere has a multiplica- 
tive effect on the target reflectivity and the convolution of two contrasting 
surface reflectances for a boundary pixel can well be approximated by their 
linear combination. So, one may define the spectral measurement of a mixed 
pixel In a wavelength band by 

Y » U Xq + (1-U) Xi (A.l) 

where Xq and Xj are the spectral measurements of two pure pixels representing 
the two clasues, say Cq and Cj, of the boundary. 

In the transformed space, discrimination between Cq and Cj is in the 
first dimension alone. Thus, It Is suffice to discuss the distribution of 
Cm for the univariate case. Suppose Xq and X\ are univariate and normally 
distributed, say Xq ~ N(-e,o^) and Xj ~ N(6,<j^). 
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It Is easy to set that the conditional distribution of Y given U • u, is normal 
with mean (l-2u)e and var1ance[u 2 + (l-u) 2 > 2 . Now writing the joint 
density, say f(y,u), of Y and U as a product of the conditional density, say 
f(y|u), and the marginal density of U, and then Integrating It with respect 
to u, the density function of Y can be expressed as follows; 


1 


f(y) -(1/20/7) 



(A.2) 


Clearly, the density function f(y) Is not of the normal type. To examine Its 
departure from normality, we next obtain Its moments and the measures of 
skewness and kurtosls. 


It can be easily verified that 
E(Y) - 0 

E(Y 2 ) - §o 2 + \ e 2 (A. 3) 

E(Y3) - 0 

E(Y 4 ) « J(7 o 4 + 8 o 2 e 2 + 0 4 ) 


The measure of skewness is zero and the measure of kurtosls Is given by 


_ 9 (7 o 4 + 8 o 2 e 2 + e 4 
5 (4 q 4 + 4o 2 "e 2 + 8 4 


To evaluate y» let 0 « o a/2, where a represents the distance between the 
distributions of classes being mixed. Then 


„ _ 9 (a 4 32a 2 + 112) 

y 5 (a 4 + iea 2 6^) 


(A.5) 
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Thus, for some typical values of A, we find y as shown below. 


A i 0 

1 

2 

3 

3.5 

4 

6 

10 

Y 1 3.15 

3.ZZ 

3.zo 

“3700“ 

“OF 

2.'75 

Z.38 

1.37 


Since y * 3 for a normal distribution, the distribution of the mixed pixel 
class Is not normal. However. Its departure from normality Is small If 
0 < a 3.5. Thus, this distribution can be approximated by normal provided 
the spectral measurements of classes constituting boundary satisfy the 
condition of a < 3.5. 
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